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The two-gluon components of the 77 and rf mesons to leading-twist accuracy 
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We critically reexamine the formalism for treating the leading-twist contributions from the two- 
gluon Fock components occurring in hard processes that involve r\ and rj mesons and establish a 
consistent set of conventions for the definition of the gluon distribution amplitude, the anomalous 
dimensions as well as the projector of a two-gluon state onto an r\ or r{ state. We calculate the 
rj, ?/ -photon transition form factor to order a s and show the cancellation of the collinear and UV 
singularities explicitly. An estimate of the lowest Gegenbauer coefficients of the gluon and quark 
distribution amplitudes is obtained from a fit to the rj, 7/ -photon transition form factor data. In 
order to elucidate the role of the two-gluon Fock component further, we analyze electroproduction 
of n, n' mesons and the g*g*n(n') vertex. 

PACS numbers: 12.38.Bx, 14.40.Aq 



INTRODUCTION 



The description of hard exclusive processes involving 
light mesons is based on the factorization of the short- 
and long-distance dynamics Q, Q- The former is rep- 
resented by process-dependent, perturbatively calculable 
parton- level subprocess amplitudes, in which the mesons 
are replaced by their valence Fock components, while the 
latter is described by process-independent meson distri- 
bution amplitudes. This work is focused on hard re- 
actions involving 77 and 77' mesons. These particles as 
other flavor neutral mesons possess SU(3)f singlet and 
octet valence Fock components and, additionally, two- 
gluon ones; to all three of them correspond distribution 
amplitudes. This feature leads, on the one hand, to the 
well-known flavor mixing which, for the 77—77' system, has 
been extensively studied (for a recent review, see 0) and, 
on the other hand, as a further complication, to mixing of 
the singlet and gluon distribution amplitudes under evo- 
lution. On the strength of more and better experimental 
data, the interest in hard reactions involving 77 and 77' 
mesons and, consequently, in the role of the two-gluon 
Fock component, has been renewed. Examples of such 
reactions are the meson-photon transition form factors, 
photo- and electroproduction of mesons or charmonium 
and £?-meson decays. 

Mixing of the singlet and gluon distribution ampli- 
tudes has been investigated in a number of papers jj]- 
[Tl|. Apart from differences in the notation and occa- 
sional misprints, different prefactors appear in the evo- 
lution kernels and in the expressions for the anomalous 
dimensions. Often the full set of conventions for kernels, 
anomalous dimensions, the gluon distribution amplitude 
and the gluon-meson projector is not provided and/or 
it is not easy to extract. This makes the comparison of 
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the various theoretical results and their applications diffi- 
cult. We therefore reexamine the treatment of the gluon 
distribution amplitude and its mixing with the singlet 
one. This analysis is performed in the context of the 777 
and 7/7 transition form factors. Applying the methods 
proposed in we calculate them to leading- twist ac- 
curacy and include next-to- leading order (NLO) pertur- 
bative QCD corrections. Our investigation enables us to 
introduce and to test the conventions for the ingredients 
of a leading-twist calculation for any hard process that 
involves 77 or 77' mesons. The most crucial test of the con- 
sistency of our set of conventions is the cancellation of the 
collinear singularities present in the parton-subprocess 
amplitude with the ultraviolet (UV) singularities appear- 
ing in the unrenormalized distribution amplitudes. Our 
analysis permits a critical appraisal of the relevant liter- 
ature H-m. 

In analogy with the analysis of the 7r7 transition form 
factor we use our leading-twist NLO results for the 
transition form factors to extract information on the 77 
and rf distribution amplitudes from fits to the experi- 
mental data [14l Il5| . In order to make contact with ex- 
periment we have to adopt an appropriate rj — 7/ mixing 
scheme. We assume particle independence of the dis- 
tribution amplitudes reducing so their number to three. 
Consequently, flavor mixing is solely encoded in the de- 
cay constants for which we use the values determined in 

m 

Our set of conventions, as abstracted from the calcu- 
lation of the transition form factor, is then appropriate 
for general use in leading-twist calculations of hard ex- 
clusive reactions involving 77 and 77' mesons. We briefly 
discuss a few of them, namely, electroproduction of the 7/ 
and 77' mesons and the vertex 3*3*77(77'), in order to learn 
more about the importance of the gluon distribution am- 
plitudes. In contrast to the transition form factors, the 
two-gluon Fock components contribute in these reactions 
to the same order of the strong coupling constant, cv s , 
as the quark- antiquark ones. The two-gluon components 
also contribute to the decays \cj — > VVjV'v'- The anal- 
ysis of these decays is however intricate since the next 
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higher Fock state of the Xcj, ccg contributes to the same 
inverse power of the relevant hard scale, the charm quark 
mass, as the cc state and has to be taken into account 
in a consistent analysis [l7|. We therefore refrain from 
analysing these decays here. 

The plan of the paper is the following: The calculation 
of the meson-photon transition form factors is presented 
in Sec.[n] In Sec. lIIII we discuss 77 — 77' flavor mixing while 
Sec. IIVI is devoted to a comparison with experiment and 
the extraction of the size of the lowest Gegenbauer coef- 
ficients of the quark and gluon distribution amplitudes. 
In Sec. [V] we investigate the role of the gluon distribu- 
tion amplitude in other hard reactions. The summary 
is presented in Sec. IVII The paper ends with three ap- 
pendices in which we compile the definitions of quark and 
gluon distribution amplitudes (App. A), calculational de- 
tails for the transition form factors (App. B) and some 
properties of the evolution kernels (App. C). 



II. THE Pj TRANSITION FORM-FACTOR 
A. The flavor-singlet case 

As the valence Fock components of the pseudoscalar 
mesons P = r\,rj ', we choose SU(3)f singlet and octet 
combinations of quark-antiquark states 1 

\qq x ) = \(uu + dd + ss) / VS) , 

\qq s ) = \{uu + dd- 2ss)/v / 6) , (2.1) 

and the two-gluon state \gg) which also possess flavor- 
singlet quantum numbers and contributes to leading 
twist. The corresponding distribution amplitudes are de- 
noted by <pp\^,g'i their formal definitions are given in 
App. ^] We emphasize that here, in this section, we do 
not make use of a flavor mixing scheme since the theoret- 
ical treatment of the transition form factors is indepen- 
dent of it. As usual the decay constants, defined by the 
vacuum-meson matrix elements of flavor-singlet or octet 
weak axial vector currents (i — 1,8) 

(0\4 5 (0)\P(p))=if P P», (2-2) 

or rather the factors f P /(2y / 2N c ), are pulled out of the 
distribution amplitudes (N c being the number of colors) . 
Hence, the quark distribution amplitudes are normalized 
to unity at any scale [i 2 

f du(f>p l (u,n 2 ) = l, (2.3) 
Jo 



1 This should not be mixed up with the usual singlet and octet 
basis frequently used for the description of r) — 77' mixing. Our 
ansatz is completely general. 



as follows from (|2.2[l and (|A9|) . From l|A10(i one has 

/ du(f>p g {u^ 2 ) = 0. (2.4) 
Jo 

There is no natural way to normalize the gluon distri- 
bution amplitude. Since the flavor-singlet quark and 
gluon distribution amplitudes mix under evolution while 
the flavor-octet one evolves independently with the hard 
scale, it is convenient to pull out of the gluon distribution 
amplitude the same factor as for the flavor-singlet quark 
one. 

As usual we parameterize the 7*(gi, /i) 7(^2, v) — * P(p) 
vertex as 

= ie 2 F Pl {Q 2 ) e^ al3 e v {q 2 ) q la q 2p , (2.5) 

where Q 2 = — q\ > is the momentum transfer, and 
Fp 1 (Q 2 ) denotes the Pj transition form factor. It can 
be represented as a sum of the flavor-octet and the flavor- 
singlet contributions 

F P ,(Q 2 ) = F^(Q 2 ) + F^(Q 2 ) , (2.6) 

where the latter one includes the quark and the gluon 
part. The leading- twist singlet contribution to order a s 
is unknown, while the octet contribution is well-known 
to this order, one only has to adapt the result for the ttj 
transitions [l^] suitably. We therefore perform a detailed 
analysis of the singlet contribution along the lines of the 
flavor-octet analysis presented in |T^ . 

For large momentum transfer Q 2 , the flavor-singlet 
contribution to the transition form factor can be rep- 
resented as a convolution (see Fig. ^ for a lowest order 
Feynman diagram) 

F^(Q 2 ) = ^=,THu,Q 2 ) ® #T(u), (2.7) 

where the symbol ® represents the usual convolution 
A(z) ® B(z) = Jg 1 dzA(z)B(z) . We employ a two- 
component vector notation 

(2.8) 

and switch to the more generic notation <f>p q = <f>p\. The 
unrenormalized quark and gluon distribution amplitudes 
<j) P r q and 4>f g are defined in Eqs. and (fXSfl . The 

parton-level subprocesses amplitudes for 7*7 — > qq, and 
7*7 ~~ > 99 are denoted by T q -q and T gg , respectively; the 
Lorentz structure is factorized out as in 1)2.5(1 . 

The distribution amplitudes §p r q and tj)^ re- 
quire renormalization which introduces mixing of 
the composite operators *S>(— z) 7+75 ft *f?(z) and 
G +a (— z) Q G a + (z). The unrenormalized distribution 
amplitude <j)p r is related to the renormalized one, <j)p, by 

(f> P r (u) = Z{u,x 1 tM 2 F )(g>(f>p{x,LJ 2 F ), (2.9) 
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FIG. 1: Lowest order Feynman diagram for the 7*7 — > P 
transition. A second diagram is obtained by interchanging 
the photon vertices. 



a matrix 2 



V qq V qg 

v gq v gg 



(2.14) 



The amplitudes T qq and T gg have well-defined expansions 
in a s , and after coupling-constant renormalization, which 
introduces the renormalization scale \j? R , they read 



« s (m|) 



4tt 



/4 



g§) T$(u) + 0(a 2 s ) 



where the UV-divergent renormalization matrix takes the 
form 



T (u) - ^ 
+ ggW - Q2 



a s {n 



4tt 



§|J T«( U ) + 0(a 2 ) 



(2.15) 



■^99 Z q g 
Zg q Zgg 



(2.10) 



The normalization factors N qq and iV sg in l|2.15|l are given 
by 



Here, /j, f represents the scale at which the singularities 
and, hence, soft and hard physics, are factorized. Owing 
to the fact that quarks and gluons are taken to be mass- 
less and onshell, T qq and T gg , calculated beyond lead- 
ing order, contain collinear singularities. The validity of 
factorization into hard and soft physics, as expressed in 
(|2.7() , requires the cancellation of these singularities with 
the UV ones from the renormalization of the distribu- 
tion amplitudes. Hence, the hard scattering amplitude 
defined by 



Tjj(x, Q 2 ,h 2 f ) = T f (u, Q 2 ) <g> Z(u, x, /4) 



(2.11) 



N qq 



N„ 



^/rifCFN qq , (2.16) 



where the flavor factor C\ takes into account the quark 
content of the qq 1 combination. It reads (see (|Al(l ) 



C x = 



el + e? + el 



(2.17) 



The number of flavors in the qq 1 is denoted by n./ and 
Cf = {N 2 — 1)/{2N C ) is the usual color factor. e a is the 
charge of quark a in units of the positron charge e. 

Inserting l|2.13[) and l|2.15|l into (|2.11|) and using (|B1|I , 
we obtain 



must be finite. Below we explicitly show this cancella- 
tion to NLO. Provided the cancellation of the singular- 
ities holds, the transition form factor can be expressed 
in terms of finite hard scattering and distribution ampli- 
tudes 



T, 



H, qq 



Q 



T (o) 

99 



47T 



Cf T nfi 



(i) / Vr 



99 \ Q2 

O(al) 



Fp 7 



(Q 2 ) 



fh 



T H {x, Q 2 ,^ ®<j> P {x,^ F ), (2.12) 



B. The NLO hard-scattering amplitude 



T, 



H,gg 



Ngg fa S (H R ) 

Q 2 



■in 



N qq 1. 



N, 



-T, 



(0) 



<."/ 



T (l) / ^R 

99 \Q 2 

V (& 

q9 Wf 



+ 0(a 2 s ) 



(2.18) 



We now proceed to the NLO calculation. The renor- 
malization matrix Z , can be shown to have the following 
form 



.(o) T (i) T (i) 
99 ' 99 ' 99 ' 
their calculation are given in App. [B] Using the results 

n(0) 

■ 99 



Results for T^> , , T gg \ and V y _and some details of 
cula 

for and V qq , it is easy to verify that 



^-yw + ^l, (2.13) 



T^{u)®V qq {u,x) = C F A ( X^)- 



(2.19) 



if dimensional regularization (D = 4 — 2e) is employed. 
Here 1 denotes the unit 2x2 matrix (with diagonal ele- 
ments S(x — u)), and the coefficient = V^(x,u) is 



2 Since we are only interested in the a s term, we suppress the label 
1 in the matrix elements of V M . 



4 



with A^Ji go being given in (|B5|) . On the other hand, where the evolution kernel V reads 



■^coi qq 1S tne res idue of the 1/e pole in Tg\ see dHJ. 

Hence, the collinear singularity present in TjQ is canceled 
by the UV singularity in Z qq and we arrive at a finite 
hard-scattering amplitude for the 7*7 — > qq subprocess 



TH,qq{x,Q 2 ,n 2 F ) 



N, 



a «(/4) r T (i) , 



(2.20) 



where 



/4 



(2.21) 



1(1) 

col,qq ' 



The quantities T^, -4^] i9f) and „4^ } are given in lB4l 



(2.22) 



defined in i|B8|) . Inserting this result into 



Next, from l|B4(l and l(B18|) . we obtain 

x col,gg( X ) ' 



lg>(«)®V„(u,«) = ^^ 



with A 



(i) 

col,gg 



(12.18(1 and taking into account (|B7() , we observe the can- 
cellation of the collinear singularity present in Tg^ with 
the UV singularity of Z qg , and we get the finite hard- 
scattering amplitude for the 7*7 — > gg subprocess 



T H ,gg(x,Q 2 ) = ^ 



where reads 



"4/4)^(1) 



47T 



T { ni g (x, Q\ Ml) = ^ §| + ^ (*) • ( 2 - 24 ) 



The functions .A^j gg and Agg ; are supplied in 1B8I 



(i) 



C. Evolution of the flavor-singlet quark and gluon 
distribution amplitudes 

We now turn to the discussion of the distribution am- 
plitude 4>p and its evolution. The matrix Z is related 
to the evolution of the distribution amplitude, and 
in (|2.13|) represents the kernel which governs the leading- 
order (LO) evolution of the flavor-singlet distribution am- 
plitude. By differentiating 1)2. 9|) with respect to \xp one 
obtains the evolution equation 0, 



V(x,u,a s (n F )) ® 4>p(u,h f ), 



(2.25) 



V = -Z 



2 d 



(2.26) 



We note in passing that the evolution equation would 
have a more complicated form if the factor fp 1 /(2y / 2N c ) 
was not pulled out of the gluon distribution amplitude. 
Inserting l|2.13[l into l|2.26|l . and using l|B2|) . one easily 
sees that 



y = y ( l) + 0{a 2 g) 

in 



(2.27) 



The results for the LO kernel are given in (|B10() 
and (|B18IB20|) . The anomalous dimensions that control 
the evolution of the distribution amplitudes can be read 
off from the relations (IC2I) : 



, </</ 



Q9 



c F 



n+1 



y/nfCF 



00 + N c 



(n+l)(n + 2) 

n(n + 3) 
3(n+l)(n + 2) 
12 

(n + l)(n + 2) 
8 

(n + l)(n + 2) 



n > 2, 
n> 2, 

n+l 



n> 2. 
(2.28) 



To leading order in a s the evolution equation (|2.25|) 
can be solved by diagonalizing the kernel V or rather the 
matrix of the anomalous dimensions. The eigenfunctions 



(2.23) can be expanded upon the Gegenbauer polynomials C n 



m/2 



with coefficients Bp^ which evolve with the eigenvalues 
7n ± ' ) of the matrix of the anomalous dimensions 

1 



~(±) = _ 
In 2 



ll q + 7 9 n 9 ± V(7^-7^) 2 + 4 7 « 



(2.29) 

The two components of the distribution amplitude <pp 
possess the expansion 



<j> Pq {x,n 2 F ) 



(j> Pg (x,n F ) 



6a: (1 — as) [l 

+ E BpM)C r ^(2x-l) 



x 2 (l-x? 



E B^ F )C^\(2x-l), 



(2.30) 

where only the terms for even n occur as a consequence 
of HA8J1 . The expansion coefficients in l|2.30|l are related 
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to those of the eigenfunctions by 



7i +) //3o 



+p^ ) 4; ) (^o) 



a a (/40 



R9 f..2j _ (+) g(+) / 2s ( a s(l4) 



+ 4n } (Mo) 



2\ \ 7i _) //3o 



(2.31) 



The coefficients Bp^ (//q) respective -Bpjf (/-to), where /Xq 
is the initial scale of the evolution, represent the non- 
perturbative input to a calculation of the transition form 
factors and are, at present, not calculable with a sufficient 
degree of accuracy. The parameters pn read 



In 



p 



(-) = t 



i 



In 



6 

7« 



In 



«(+) ' 

(2.32) 

We note that the anomalous dimensions satisfy the re- 
lation 



^99 
/n 



7ri 



(±) 



/n 



.</.</ 



,(±) 



mi 



(2.33) 



in 



In - 7, 

Comparison of (|2.28|) and (|2.29jl reveals that 7^ 

for all n and 7„ — > 7^ for n — * 00. 

It is important to realize that any change of the defini- 
tion of the gluon distribution amplitude <|A5I) is accom- 
panied by a corresponding change in the hard scattering 
amplitude. Suppose we change 4>p g by a factor a 



,J Pg 



P g 



(2.34) 



Since any physical quantity, as for instance the transi- 
tion form factor, must be independent of the choice of 
the convention, the projection (|A14|) of gg state onto a 
pseudoscalar meson state is to be modified by a factor 
1/ct, i.e., 



-pger 



I 



V 9 



(2.35) 



and the hard-scattering amplitude becomes altered ac- 
cordingly. As an inspection of Eqs. I|2.30[l - i|2.32[l reveals, 
the change of the definition of the gluon distribution am- 
plitude l|2.34Jl has to be converted into a change of the 
off-diagonal anomalous dimensions and the Gegenbauer 
coefficients Bp^ in order to leave the quark distribution 
amplitude as it is: 



qg,a _ _ qg 
In In ' 

a 



In w In > 



(2.36) 



and 



B 



( - )CT (^)=a4-)(^), 



B 



(+)" 



04) 



B P + J (Mo) 



(2.37) 



implying 



'Pn KH'FJ — a B Pn (pF) 1 LJ l'n >/'/.'' '->!•>, 

(2.38) 

We finally mention that, as can be easily seen from Eq. 
(12.34(1 and the evolution equation H2.25J1 . along with the 
change of the anomalous dimensions (|2.3(i|) the kernels 



I), 



V qg and V gq become modified. 

The results for the anomalous dimensions can also be 
understood in the operator language, i.e., by consider- 
ing the impact of a change of the definition of the glu- 
onic composite operator on the anomalous dimensions 
(for comments on the use of the operator product expan- 
sion, see, for instance, Refs. pQ)- One finds that only 
the anomalous dimensions 7^ ff and 7,^ become modified, 
while the diagonal ones and the product 7 r ? 9 7^ 9 , and 
consequently the eigenvalues 'Jn 1 remain unchanged. 
Redefinition of the gluonic composite operator implies 
a corresponding change of the gluon distribution ampli- 
tude. 

We are now in the position to compare the results pre- 
sented in this work with other calculations to be found 
in the literature. The entire set of conventions is not al- 
ways easy to extract from the literature since often only 
certain aspects of the flavor-singlet system are discussed. 
For instance, in Ref. [9j only the evolution kernels are 
investigated, or in Ref. Q only the anomalous dimen- 
sions. Using results from such work in a calculation of a 
hard process necessitates the use of corresponding con- 
ventions for the other quantities. Care is also required 
if elsewhere determined numerical results for the Gegen- 
bauer coefficients B p J or Bp n are employed since, ac- 
cording to (|2.37() and l|2.38[l . they are convention depen- 
dent. For future reference, we systematize in Tab. [I] the 
important ingredients for the three conventions encoun- 
tered in the literature. Our expressions for the kernels 
and the anomalous dimensions correspond to the ones 
obtained in Q (up to a typo in V gg ). In Refs. [H EH 
the anomalous dimensions controlling the evolution of the 
forward and non-forward parton distribution were stud- 
ied to NLO. Since the non-diagonal anomalous dimen- 
sion for the odd parity case coincides with our ones |l9j| , 
we observe that the convention a — *J CfJti} is used in 
|ToL ITlT | . The only result we do not understand is the one 
presented in Ref. pg: There is an extra factor of 1/2 in 
V gq which changes the product of prefactors. Moreover, 
there are factors 1/3 and 3 apparently missing in j qg and 
-f gq . We note that occasionally the factor [x(l — x)]" 1 
appearing in our projector (| A 1411 is absorbed into the 
gluon distribution amplitude 0, Q • This arrangement is 
accompanied by corresponding changes of the evolution 
kernels, see QB21jl . 

Although, from the point of view of derivation, the 
conventions which lead to l|2.28[) and l|A14() seem to be 
the most natural ones, it is perhaps more expedient to 
use the same conventions for the anomalous dimensions 
as for polarized deep inelastic lepton-proton scattering 
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TABLE I: List of common conventions for the anomalous di- 
mensions and the gg projector. Quoted are the prefactors of 
the non-diagonal anomalous dimensions 112.2811 and of the gg 
projector 1 A 1411 for various choices of a in 12.351 12.361 . We 
also list references where these conventions for the anomalous 
dimensions are used. 



In 


^31," 
In 


■pga 


references 










\Jnf Cf 


sjn s Cf 


1 


1 


C f 


n f 


fc£ 








V n f 


n f 


C F 


V Cf 


[9.UJ 



[20|. which correspond to 




(2.39) 



The corresponding set of conventions will be used in the 
rest of the paper. The non-diagonal anomalous dimen- 
sions then read 



^13 

In 



In 



C F 
n f - 



n(n + 3) 



3(n + l)(n + 2) 
12 



(n+l)(n + 2) 
and the gluonic projector 

Sab 



n>2, 
n>2, (2.40) 




Cf 

n f ^N 2 - 1 u(l-tt) 



(2.41) 



Along with these definitions, Eqs. (|2.30() - (|2.32|l have to 
be used. 

To the order we are working, the NLO evolution of the 
quark distribution amplitudes should in principle be in- 
cluded (the convolution of the NLO term for <fip g with 
Tn,gg contributes to order a 2 s ). To NLO accuracy the 

3/2 

Gegenbauer polynomials C n are no longer eigenfunc- 
tions of the evolution kernel, so that their coefficients 
B Pn do not evolve independently 0, |2l|. In analogy 
with the pion case (2^ , the impact of the NLO evolution 
on the transition form factors is expected to be small 
compared with the NLO corrections to the subprocess 
amplitudes. Therefore we refrain from considering NLO 
evolution. 



D. The NLO result for the transition form factor 

To end this section we quote our final result for the 
flavor-singlet contribution to the P7 transition form fac- 
tor to leading- twist accuracy and NLO in a s . The re- 
sult, obtained by inserting H2.20|) and l|2.23|l (multiplied 



by a 1 = \/ Cf /n / according to the new normalization 
of the gluonic projector) into (|2.12|) . is 



fpc 



^ {^-(aO®<M*,/4) 



Air 



(2.42) 



A subtlety has to be mentioned. The singlet decay con- 
stant, fp, depends on the scale but the anomalous di- 
mension controlling it is of order a 2 23]. In our NLO 
calculation this effect is tiny and is to be neglected as 
the NLO evolution of the distribution amplitude. 

For completeness and for later use we also quote the 
result for the flavor-octet contribution to the P7 transi- 
tion form factor at the same level of theoretical accuracy. 
In our notation it reads 



fpCs 



C F TgJx, Q 2 , &) ® 0p 8 (x, A) 



4tt 



(2.43) 



where the renormalized hard scattering amplitude is 
given in (|2.2Q(I and the charge factor Cs is obtained with 
the help of (|Alf) 



C H = 



2e 



V6 



(2.44) 



The octet distribution amplitude, 4>ps, being fully anal- 
ogous to the pion case, has the expansion 



<f>P S (x, fj.%) = 6x(l- x) [1 



n=2A.. 



BZ n {p? F )C*/\2x-l) 



(2.45) 

where the Gegenbauer coefficients evolve according to pj 



7„"//3o 



(2.46) 



Summing the flavor-singlet and octet contributions ac- 
cording to l|2.6f) . we arrive at the full transition form fac- 
tors for the physical mesons. 

As has been pointed in Refs. 0, IT3L |24| . in the limit 
Q 2 — > 00 where the quark distribution amplitudes evolve 
into the asymptotic form 



(f>As(x) = 6x(l - x) 



(2.47) 
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and the gluon one to zero, the transition form factor be- 
comes 



F 



Pi 



Q 2 



Q 2 



3 7T 



(2.48) 



/J> ff combines the decay constants with the charge factors 



rcS 
JP 



1 

7! 



fp + 2V2/J, 



(2.49) 



The result (|2.48|) holds also for the case of the pion with 
f° s replaced by / w . In Q an interesting observation has 
been reported: if the transition form factors for the ir, 77 
and 77' are scaled by their res pec tive asymptotic results, 
the data for these processes 0, fall on top of each 
other within experimental errors. This can be regarded 
as a hint at rather similar forms of the quark distribution 
amplitudes in the three cases and a not excessively large 
gluon one. 



III. 77-77' MIXING 

Using the results l|2.42|l and (|2.43l) for the transition 
form factors, one may analyze the experimental data ob- 
tained by CLEO [lj and L3 [H with the aim of ex- 
tracting information on the six distribution amplitudes 
4>pi{x, /!q), i = 1,8, g or rather on their lowest Gegen- 
bauer coefficients Sp n (/ig)- In principle, this is an ex- 
tremely interesting program since it would allow for an 
investigation of 77 — 77' flavor mixing at the level of the 
distribution amplitudes. In practice, however, this pro- 
gram is to ambitious since the present quality of the 
data is insufficient to fix a minimum number of six coef- 
ficients which occur if the Gegenbauer series is truncated 
at 77 = 2. Thus, we are forced to change the strategy and 
to employ a flavor mixing scheme right from the begin- 
ning in order to reduce the number of free parameters. 

Since in hard processes only small spatial quark- 
antiquark separations are of relevance, it is sufficiently 
suggestive to embed the particle dependence and the mix- 
ing behaviour of the valence Fock components solely into 
the decay constants, which play the role of wave func- 
tions at the origin. Hence, following 0, [24|, we take 

4>Pi = <$>i ! (3-1) 

for i = 8,1, g. This assumption is further supported by 
the observation [24],|25| that, as for the case of the pion 
[iH. I22I |26| , the quark distribution amplitudes for the 77 
and 77' mesons seem to be close to the asymptotic form 
4>as(x) for which the particle independence (|3.1|) holds 
trivially . Note that we switch now back to the origi- 
nal notation for the singlet distribution amplitude intro- 
duced in Sec. Ill Al 



9 Pi 



?Pq ■ 



B Pn 



D Pr, 



(3.2) 



The decay constants can be parameterized as [16|, |2 



f* = f 8 COSt 



ft = / 8 sin( 



fj, =-/isin0 1: 
fh = fi cos 0i 



(3-3) 



Numerical values for the mixing parameters have been 
determined on the basis of the quark-flavor mixing 
scheme Ha : 



h 
h 



1.26U 
1.17/* 



-21.2° 
-9.2° . 



(3.4) 



The value of the pion decay constant is / w = 0.131 GeV. 
As observed in (see also Q) f] — ff flavor mixing can 
be parameterized in the simplest way in the quark-flavor 
basis. The mixing behaviour of the decay constants in 
that basis follows the pattern of state mixing, i.e. there 
is only one mixing angle. The basis states of the quark- 
flavor mixing scheme are defined by 



\r) q ) = cos tp\q) + sin tp\rf), 
I ?7 S ) = ~smip\i]} + cos ip\r]'}, 



(3.5) 



and the strange and non-strange decay constants are as- 
sumed to mix as 



1 



J ' q cos (p, /* = -f s sinip, 

f q sin ip , = f a cos Lp . (3.6) 

As demonstrated in |16| this ansatz is well in agreement 
with experiment. The occurrence of only one mixing an- 
gle in this scheme is a consequence of the smallness of OZI 
rule violations which amount to only a few percent and 
can safely be neglected in most cases. SU (3) f symmetry, 
on the other hand, is broken at the level of 10 — 20% as 
can be seen, for instance, from the values of the decay 
constants fs and f±, and cannot be ignored. 

Using p.lfl and particle independence, we obtain for 
the valence Fock components of the basis states l|3.5ll 

h 



Wq) = 



.Tyy- W x >/4') \qq) + 0o PP (x,hf) \ss) 

^2/3 <t> g (x,n F ) \gg) 



f S 



[<l>opp(x, Vf) 111) + fafaA) 



+<t> g {x,n F ) \gg)/Vs 



(3.7) 



where qq is short for the combination {uu-\-dd)/y2 and 

4> q = \{<h + 2( M i (t>s = I (2<^8 + 4>i) > 



V2, 
3 



(3.8) 



In deriving (|3.7J) we made use of the relations 

1 /, , n . 1 f 



cos(ip - 6 S ) 

sin(<p-0 8 ) = l/rTi sin((/3-0i) 



2/- _ S _/.. ' 2 /■/ 



3/8 



3/1 ' 

(3.9) 
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which can readily be obtained from results on decay con- 
stants and mixing angles reported in |16|. 

In (|3.5|) the ss (qq) Fock component appears in the r\ q 
(j] s ). These respective opposite Fock components lead to 
violations of the OZI rule if they were not suppressed. In 
order to achieve the mixing behaviour 13. 5[) . (|3.6[) and, 
hence, strict validity of the OZI rule, 
which implies 



^opp 



must be zero 



4>&{x,Hf) = 4>i{x,fi%) = (f> q (x,(j, 2 F ) = (f>s{x,n 2 F ) . (3.10) 

However, except the distribution amplitudes assume the 
asymptotic form, this can only hold approximately for a 
limited range of the factorization scale since the evolution 
of the distribution amplitudes will generate differences 
between <f>i and <ps and, hence, the respective opposite 
Fock components. In order to guarantee at least the ap- 
proximate validity of the OZI rule and the quark-flavor 
mixing scheme as is required by phenomenology, we de- 
mand in our analysis of the transition form factor data 
that 



<I>as{x) 



« 1, 



(3.11) 



for any value of x. 



IV. DETERMINATION OF THE 
DISTRIBUTION AMPLITUDES 

Before we turn to the analysis of the P7 transition 
form factor data 0, 0] and the determination of the 
r\ and 77' distribution amplitudes a few comments on the 
choice of the factorization and renormalization scales are 
in order. A convenient choice of the factorization scale 
3 is n 2 F = Q 2 , it avoids the \n)i 2 F /Q 2 terms in (12 .20(1 
and (|2.23|) . Another popular choice is fj, F = Q 2 /2 which 
reflects the mean virtuality of the exchanged quark. This 
choice facilitates comparison with the pion distribution 
amplitude as determined in 0] in exactly the same way 
we are going to fix the rj and rf distribution amplitudes. 
For the renormalization scale we choose n 2 R = Q 2 /2 for 
which choice arguments have been given on the basis of 
a next-next-to-leading order calculation of the pion form 
factor : 12j. 

The transition form factor is evaluated using the two- 



MS 



loop expression for a s with four flavors and A 
305 MeV |28j. The numerical values for the decay con- 
stants and mixing angles arc given in 13. 4|) . As the start- 
ing scale of the evolution we take /z 2 , = 1 GeV 2 . 

A comparison of the leading-twist NLO results evalu- 
ated from the asymptotic quark distribution amplitudes 



3 A detailed discussion of the the role of the factorization scale 
and the resummation of corresponding logs is presented in Refs. 



> 

o 



o 



o 



0.3 



0.2 



0.1 




B!=-0.04 



,(LO) 
asym. DA (NLO) 

B' = 21 B^-0.04 

B' = -3 B' =-0.12 



10 
Q 2 



15 



20 



FIG. 2: The scaled Pj transition form factor vs. Q 2 . Dotted 
(long-short dashed) lines represent the LO (NLO) predictions 
for the asymptotic distribution amplitudes. Solid (dashed) 
lines are results obtained with Bf (/.to) = 21 (-3), BgO-'o) = 
-0.04 (-0.12) and Bf(^) = -0.04 (fj, 2 F = Q 2 , fi% = Q 2 /2, 
(Xq = 1 GeV 2 ). The shaded areas indicate the range of the 
NLO predictions for B 2 and Bf inside the allowed region (see 
text). Data taken from (THHEl (rhombs represent the Q 2 F, ] i 1 
data, squares the Q 2 F n ~, ones). 



(|2.47|l (the gluon distribution amplitude is zero in this 
case) with experiment [Til Il5| is made in Fig. [21 It 
reveals that the distribution amplitudes cannot assume 
their asymptotic forms for scales of the order of a few 
GeV 2 ; the prediction for the case of rj' lies about 10% 
above the data. This parallels observations made for the 
case of the n^f transitions [T3L |22| . 

Next let us inspect the Gegenbauer expansion of the 
transition form factor. For x-independent factorization 
and renormalization scales the integrations involved in 
(|2.42l) and l|2.43|) can be performed analytically leading 
to the expansion 



F Pl (Q 2 ) 



6/>Ci 



1 



B\{n F 



5 a s (/4) 

3 7T 



10487 



59 
72 



5, Q 2 

c lll_ 2" 

6 n F 



4500 
55 

1296 ~ 



+M 04) ( 

4 KH'FJ I 10125 



91 Q 2 

m -^r- 

75 n 2 F 

J_i 

108 ix 2 F 

_7_, ^ 

675 n n 2 F 



+ 



(4.1) 



Particle independence of the distribution amplitudes is 
used in this expansion. A similar expansion holds for the 
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octet contribution with the obvious replacements fj, — > 
ff,, B 1 B% and B£ -> 0. The expansion of the octet 
contribution is analogous to that one of the transition 
form factor [H[l3|. 

In the expansion (|4.1|l one notes a strong linear corre- 
lation between B\ and B\, only the mild logarithmic Q 2 
dependence due to evolution and the running of a s re- 
stricts their values to a finite region in parameter space. 
The gluon contributions to the form factors are strongly 
suppressed, they appear only to NLO and the numeri- 
cal factors multiplying their Gegenbauer coefficients are 
small. The coefficients flf an d B\ are also correlated. 

With regard to these correlations and in view of the 
errors of the experimental data [T3. ITfij as well as the 
rather restricted range of momentum transfer in which 
they are available, we are forced to truncate the Gegen- 
bauer series at n = 2. Truncating at n — 4 does not lead 
to reliable results in contrast to the simpler case of the 
pion where this is possible [l3j. A fit to the CLEO and 
L3 data for Q 2 larger then 2 GeV 2 provides 



#f0*o) = -0.04 ±0.04, 
B l(l4) = -0.08 ±0.04, 
BZ(j£) = 9 ±12, 



(4.2) 



where the values of the Gegenbauer coefficients are ob- 
tained for the factorization scale fip = Q 2 . We repeat 
that /j,q = 1 GeV 2 and the gluonic Gegenbauer coefficient 
is quoted for the normalization a — \Jn$ /Cf- For com- 
parison we also determine the Gegenbauer coefficients for 
fip = Q 2 /2; the values found agree with those quoted in 
(|4.2|) almost perfectly. The quality of the fit is shown in 
Fig. |2 The coefficients B\ and £?f are strongly corre- 
lated as can be seen from Fig. |3| The results H4.2(l satisfy 
VUl-Bf (/40 - #2 04)1 /3 < 0-02 for all (i% > /z§. This 
meets the requirement and, therefore no substan- 

tial violations of the OZI rule follow from our distribution 
amplitudes. It moreover implies the approximative va- 
lidity of the quark-flavor mixing scheme advocated for in 
Ref. ^(J- I n Fig- 01 we present the singlet and gluon dis- 
tribution amplitudes at the scale /z 2 , obtained using the 
face values from 1)4.2(1 . Both amplitudes are end-point 
suppressed as compared to the asymptotic one. This 
property holds for all values of B\ and f?f inside the 
allowed region (|4.2|l . 

The values of B\ and B\ agree with each other within 
errors as well as with the Gegenbauer coefficient (/tig) 
of the pion distribution amplitude for which a value of 
— 0.06±0.03 has been found in |l3| from an analysis along 
the same lines as our one. Thus, the three quark distri- 
bution amplitudes are very similar. This result explains 
the observation made in |3( and mentioned by us at the 
end of Sec. II D that the data on three transition form 
factors fall on top of each other within errors if the form 
factors are scaled by their respective asymptotic results 
(|2.48[l . The rjcj transition form factor, on the other hand, 
behaves differently |2^. The rj c mass provides a second 



30 



20 



10 
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-0.15 



-0.10 



-0.05 



Bl 



FIG. 3: la x 2 -contour plot for the coefficients B\{iM)) and 
B% (Mo) obtained from a three-parameter fit to the CLEO and 
L3 data on the r/, »/~7 transition form factors. Values of the 
Gegenbauer coefficients refer to Hq = 1 GeV 2 ; the factoriza- 
tion scale is /x% — Q 2 . 




x, (4), B^-0.08 
x, £), B== 9 



0.5 

x 



FIG. 4: Flavour-singlet and gluon distribution amplitudes at 
the scale /xq = 1 GeV 2 obtained using the face values B\ 
and Bf from 114.211 . The asymptotic distribution amplitude is 
included for comparison. 



large scale which cannot be ignored in the analysis |30|. 

We emphasize that our results on the r\ and T]' dis- 
tribution amplitudes are to be considered as estimates 
performed with the purpose of getting an idea about the 
magnitude of the gluon distribution amplitude. As has 
been discussed in detail for the case of the ttj transition 
form factor in [l3j, allowance of higher Gegenbauer co- 
efficients in the analysis will change the result on B^ , 
essentially the sum of the B% is fixed by the data on the 
transition form factor. This ambiguity also holds for the 
case of the rj and 77'. Taking a lower renormalization scale 
than we do which may go along with a prescription for 
the saturation of a s and thus including effects beyond 
a leading-twist analysis, will also change the results for 
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A 



\k' = k + A 





msmwrn 



p' =p + A 



FIG. 6: Sample leading order Feynman diagrams that con- 
tribute to the subprocess amplitude j^q — > qq i q. 



FIG. 5: The handbag- type diagram for meson electroproduc- 
tion off protons. The large blob represents a generalized par- 
ton distribution, the small one the subprocess ^* L q — > P q. 
The momentum transfer is t = A 2 . 



the Gegenbauer coefficients. Another source of theoret- 
ical uncertainties in our analysis is the neglect of power 
and/or higher-twist corrections. Thus, for instance, in 
Refs. j2J,|2£j the LO modified perturbative approach 
has been applied where quark transverse degrees of free- 
dom and Sudakov suppressions are taken into account. 
In this case the asymptotic distribution amplitudes lead 
to good agreement with the data on the transition form 
factors. 



V. COMMENTS ON OTHER HARD 
REACTIONS 

In this section we make use of the results obtained in 
the preceding sections and calculate other hard processes 
involving rj and rj mesons in order to examine the role of 
the gg Fock component further. 



A. Electroproduction of 77,7/ mesons 

As a first application of the gluon distribution ampli- 
tude extracted from the 777 and 7/7 transition form fac- 
tors we calculate deeply virtual electroproduction of jj 
and rj mesons off protons. It has been shown j33, Hij 
that for large virtualities of the exchanged photon, Q 2 , 
and small momentum transfer from the initial to the fi- 
nal proton, t, electroproduction of pseudoscalar mesons 
is dominated by longitudinally polarized virtual photons 
and the process amplitude factorizes into a parton-level 
subprocess 7^ q — > Pq and soft proton matrix elements 
which represent generalized parton distributions |34| , see 
Fig. [SJ The meson is generated by a leading-twist mech- 
anism, i.e., by the transition qq — > P mediated through 
the exchange of a hard gluon. For the production of 77 
and 77' mesons, however, one has to consider the gluon 
Fock component as well which, in contrast to the case of 
the transition form factors, contributes to the same order 
of a s as the qqi components. The gluonic contribution 
has not been considered in previous calculations of the 
electroproduction cross sections [35l . 

The helicity amplitude for the process "f^p — > Pp is 



again decomposed into flavor octet and singlet compo- 
nents, qqi — > P 



M%1 0± = ^ee.Qv/l - ^ 



dx 



■ H 



Pi 

o±,o± 



, (5.1) 



where H a and E a are the axial vector and pseudoscalar 
generalized parton distributions for the emission and re- 
absorption of quarks of flavor a. The C* are flavor fac- 
tors for the qqi components of the meson P; they can 
be read off from I|A1|1 . The quark subprocess amplitudes 
7^ q — ► qqi q are calculated from the LO Feynman dia- 
grams for which examples are shown in Fig. [H] 36] 



QV~- 



Q 2 



< I dr 
'0 



(1 - r)Q 2 - rt 
u It' 

S 1 — T U 



(5.2) 



They are expressed in terms of the subprocess Mandel- 
stam variables s,u,t = t where s + t + u = —Q 2 , and 
hold for any value of Q 2 and t. For the deeply virtual 
kinematical region of large Q 2 and — t <C Q 2 , it is more 
appropriate to use the scaling variables £ and x. The 
skewness is defined by the ratio of light-cone plus com- 
ponents of the incoming (p) and outgoing (p') proton 
momenta 



(p-p') 
{p + p')+ 



(5.3) 



For large Q 2 the skewness is related to x-Bjorken by £ ~ 
XBj/2- The average momentum fraction the emitted and 
reabsorbed partons carry, is defined as 



(fc + fcQ- 
(p + p')- 



(5.4) 



Here, k and k' are the momenta of the emitted and re- 
absorbed partons, respectively. For —t Q 2 the Man- 
delstam variables are related to the skewness and the 
average momentum fraction 
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2£ 



(2 - , 



2i 



(5 + 



(5.5) 



Rewriting the subprocess amplitude in terms of £ and 
x and inserting the result into the factorization formula 
(|5.1|l , one arrives at the well-known result for the leading- 
twist contribution to deeply virtual electroproduction of 
pseudoscalar mesons 



47ra s (/x|j) C F 



± 



/p^I 3 ^ I dr 
Jo 



<t>i(T, fi%) 



dx 

e 



l£ 



(5.6) 



Next we calculate the subprocess amplitude for the 
gluonic component of the meson, q — ► gg q. There 
are six graphs that contribute to the subprocess. Three 
representative ones are depicted in Fig.Q the other three 
ones are obtained from these by interchanging the gluons. 
We find for that subprocess amplitude the result 



^0±,o±(^! Q 2 ) 



_j_ A f 2 \ fh CF Q 



-t 



1 dT M T > v-p) 



T 2 (l 



+ 8 y/—US 

(5.7) 



In deriving this expression we made use of the antisym- 
metry of the gluon distribution amplitude l|A8(l . The glu- 
onic contribution to the ^* L p — > Pp helicity amplitudes 
reads 



M?l 0± = 



dx 



Pg 

o±,o± 



(5.8) 



The full 7£ p — > Pp amplitudes are the sum of the flavor 
octet and singlet contributions 1|5.6[) and the gluonic one 
(|5.8() . In the deeply virtual region, however, the gluon 
contribution is suppressed by t/Q 2 as one readily ob- 
serves from (|5.7f) . It is, therefore, to be considered as a 
power correction to the leading quark contribution (|5.6|) 



and is to be neglected in a leading-twist analysis of deeply 
virtual electroproduction of rj and rj' mesons. 

One may also consider wide-angle photo- and electro- 
production of r\ and rj' mesons. Using the methods pro- 
posed in |37| for wide-angle Compton scattering, one can 
show that for wide-angle photo- and electroproduction 
of pseudoscalar mesons the factorization formulas i|5.1|) 
and l|5.8l) hold as well provided —t and — u are large as 
compared to the square of the proton mass and Q 2 -C —t 
|36| . To show that one has to work in a symmetric frame 
in which the skewness is zero. One can also show that, in 
this situation, s and u are approximate equal to the Man- 
delstam variables for the full process, s and u, respec- 
tively. Thus, in the wide-angle region and for Q 2 <C —t, s 
but non-zero, (|5.1() and (|5.8() simplify to 

M^ Q± (s 1 t 1 Q 2 «-t) = en^ 0± J2^C n a R a A (t), 

a 

M P | i0± ( S ,t,Q 2 «-<) = e7^ i0± ^e o i^(t), 

a 

(5.9) 

where the form factors R a A represent 1/x moments of 
the generalized parton distributions H a at zero skew- 
ness. These form factors also contribute to wide-angle 
Compton scattering |37j |. The amplitudes for transver- 
sally polarized photons can be obtained analogously. In 
contrast to the case of deeply virtual electroproduction 
[38| , factorization for these amplitudes holds in the wide- 
angle region, too. 

In order to estimate the size of the gluon contribution 
to wide-angle electroproduction of 77, 77' mesons, we plot 
in Fig. |H1 the ratio 



M 



Pg 

o±,o± 



-r 



M P1 



o±,o± 



2s 2 + t 2 + ts 



dr 



<t>g{T, Mf) 

r 2 (l-r) 



1 dr Mr, Mf) 



(5.10) 

evaluated from the distribution amplitudes (|4.2|l for 
which the ratio of the integrals is ~ —5 S|(/x 2 7 )/18. The 
ratio may be large in particular in the backward hemi- 
sphere. Thus, at least for electroproduction of rj mesons 
the gg Fock component should be taken into account 
for sufficiently large momentum transfer. For the pro- 
duction of the r\ meson it plays a minor role since r\ 
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FIG. 8: Ratio of gluon and flavor-singlet quark amplitudes 
for wide-angle electroproduction of r\ or vj mesons (/if. = 
10 GeV 2 ). The shaded area indicates the range of predictions 
evaluated from B\{^) and Bf(/io) inside the allowed region 
according to Fig. |3] 



a) 



b) 



^ 1 



FIG. 9: Relevant lowest order Feynman diagrams for the 
9*9* —> <M ( a ) and 9*9* — * 55 subprocess (b). 



production is dominated by the flavor-octet contribution 
(/rj/ ' fri = 0.16). Note, however, that the normalization 
of the meson electroproduction in both the regions, the 
deeply virtual and the wide-angle one, is not well under- 
stood in the kinematical region accessible to present day 
experiments. 



B. The g*g* P vertex 

A reliable determination of the g*g*rj' vertex is of im- 
portance for the calculation of a number of decay pro- 
cesses such as B — ► rj'K, B — > v'N s , or of the hadronic 
production process pp — > rj X. The q*q*ri' vertex has 
been calculated by two groups recently [3!j 0] ■ We re- 
analyze this vertex to leading-twist order using our set of 
conventions. This will allow us to examine the previous 
calculations, and provide predictions for Pg* transition 
form factor using the Gegenbauer coefficients 14.21) in the 
distribution amplitudes. 

We define the gluonic vertex in analogy to the electro- 
magnetic one , see (|2.5|l . as 



The calculation of the transition form factor to leading 
twist accuracy and lowest order in a s parallels that of the 
meson-photon transition form factor which we presented 
in some detail in Sec. |n] In contrast to the electromag- 
netic case, however, already to the lowest order in a s 
the two partonic subprocesses g*g* — > qq and g*g* — > gg 
contribute. The relevant Feynman diagrams are shown in 
Fig. 1^1 There are a few more diagrams which involve the 
triple and quadruple gluon vertices. The contributions 
from these diagrams are separately zero when contracted 
with either the qq or the gg projectors (|A1 Lf) . H2.41fl . 
The following result for the Pg* transition form factor 
can readily be obtained 



_> \ fp 



F Pg ,(Q ,u) = 47ra s (/4) ^ 



Nr. 



(5.13) 



where 



1 ab 



fiisot/3 



(5.11) 



where q\ and q 2 denote the momenta of the gluons now 
and a and b label the color of the gluon. It is evident 
that the transition to a colorless meson requires the same 
color of both the gluons. We consider space-like gluon 
virtualities for simplicity; the generalization to the case 
of time-like gluons is straightforward. We introduce an 
average virtuality and an asymmetry parameter by 



Q 2 = -\{<ll + fa- 



ll 



if 



4 



(5.12) 



The values of u> range from —1 to 1, but due to Bose 

symmetry the transition form factor is symmetric in this 

2 2 

variable: F Pg *(Q ,oj) = F Pg ,(Q ,-v). 



dx 



<t>g{x,H? F ) 



1 - 2x 



XX 



1-oj 2 {\-2x) 2 



(5.14) 



There is no contribution from the qq s component to this 
vertex. 

Inserting the Gegenbauer expansions i|2.3U|) into (|5.14|) 
the integrals can be performed analytically term by term 
analogously to l|4.1|) resulting in the expansions 



Aqq{L0) = Cq(uj) + C 2 (uj) #2 (Mf) + ■ 

A gg (uj) = g 2 (w) Bf(/4) + — , 



(5.15) 
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FIG. 10: Functions Co, C2, and (?2, denned in Eq. I5.16L 



where 
c (lu) 

C2(w) 

92 (w) 



3 
3 

-5 
I2u 4 



1 , 9N , 1 + LU 

1 fl-w ) n- 

2lo V ; 1 



15 - 13w 2 (5 - 6lu 2 + lu 4 ) In 

2lo v ' 1 — w 



2w 2 



3 

2^ 



(l-c 2 )ln T 



1 + w 



(5.16) 



The behaviour of functions Cq(u>), C2(lo), and <?2 (u>) is 
illustrated in Fig. 1101 Examining the function C2(u>) and 
Eq. Q5.15)) . one notice that the form factors become in- 
creasingly less sensitive to the coefficients B\{^) with 
decreasing \lo\. This behaviour is characteristic of all 
functions c n (u) (n > 0) |13j. On the other hand, the 
functions cq{lu) and 52(1^) do not depend so drastically 
on lu and they are non-zero at w = 0. One can easily 
show that all g n (cu), for n > and even, possess this 
property. 

Let us discuss two interesting limiting cases. For 10 <C 
1, i.e., for qf w q 2 , the form factors behave as 



F Pg *(Q' 



V3Q 2 



fh 



1 

12 



0{lu 



4 a 2 ) 



(5.17) 



Thus, the limiting value for lu — > is sensitive to the 
form of the gluon distribution amplitude while it does 
not depend on the Gegenbauer coefficients of the quark 
one. This is to be contrasted with the P7* transition 
form factor which, according to |13| . is independent of 
both the quark and the gluonic Gegenbauer coefficients 
in the limit lu — *■ 0. 

For lu — > ±1, i.e., in the limit where one of the gluons 



goes on-shell, the Pg transition form factor becomes 
Ay/Zna s {n 2 R )) 



F Pg (Q 2 ,Lu 



±1 



fh 



(5.18) 



where Q 2 = —q 2 (—q 2 ) as in the electromagnetic case. In 
Fig.^Jwe display our predictions for the scaled r/'g* tran- 
sition form factor evaluated from the distribution ampli- 
tudes determined in Sec. II VI choosing /uf* = n\ = Q . 
Given the large difference in the magnitude of B\(ijfy 
and B°(p 2 ), see ijOl . we observe a strong sensitivity 
of the Pg transition form factors on the gluon distribu- 
tion amplitude in contrast to the electromagnetic case. 
Due to the badly determined coefficient B^ the uncer- 
tainties in the predictions for F^g* are large. Because 
of the smallness of the mixing angle 9i, see (|3.3ll and 
(|3.4|) . the rig* transition form factor is much smaller 

then the rj'g* one. The ratio of the two form factors, 

2 2 

Fr, n *(Q , lu) I F^i fl . (Q ,lu) is given by — tan(V This re 



suit offers a way to measure the angle 0\ as has been 
pointed out in 16] . 

Let us compare our results for the nj7* transition form 
factors with those presented in Refs. [3ElEi. First we re- 
mark that there is perfect agreement for the contribution 
from the meson's qq 1 component. As for the contribution 
from the gluonic component we differ by a factor l/(2n/) 
from Refs. [3!jE(j 4 . Furthermore, in (43, there is an ad- 
ditional factor of lu multiplying the gluonic term render- 
ing it antisymmetric in lu in conflict with Bose symmetry. 
We suspect that a gluonic projector ~ e^ va P qi a q2p / Q 
is used in |4Cj which turns into ~ lue^ in a frame where 
the meson moves along the 3-direction. This is in conflict 
with HA12|) . l|A13jl except at lo = 1. 

The origin of the missing factor l/(2n/) is not easy to 
discover since in Refs. 39, 40] the form of the gluonic pro- 
jector is not s peci fied. Given the anomalous dimensions 
quoted in |39l l40j|. which are the same as in l|2.40|l . this 
incriminated factor cannot be assigned to a particular 
normalization of the gluonic projector, (|2.41l) must be ap- 
plied. On the other hand, using a = 1/(2-^/71/ Cp) as the 
normalization of the gluonic pro ject or, the results for the 
transition form factors given in |39t |4(J would be correct 
(ignoring the problem with the factor lu in |40j), provided 
the corresponding anomalous dimensions are applied, see 
(|2.36l) . and they differ from the ones quoted in these pa- 
pers. Hence, the quoted anomalous dimensions and the 



4 We corrected a typo in l39l where only the case of ui = 1 has been 
dealt with - the relative sign between the contributions from the 
two Feynman diagrams shown in Fig. |5] (b) should be minus. 
Moreover, in this work Ohrndorf's results Q for the anomalous 
dimensions are used which are flawed while they have the same 
normalization as in I2.4UI . 
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FIG. 11: Predictions for the rf ' g* transition form factor as a 

2 

function of lo for two values of Q . The shaded areas indicate 
the range of predictions evaluated from -82(^0) an< l ^Kmo) 
inside the allowed region according to Fig. 



result for the gluon part of the hard-scattering amplitude 
seem not to be in agreement. 

In Ref. the leading term of the expansion (I5.17|) 
has been derived from the results presented in |4(J] and 
it therefore disagrees with our result. 



As an application of the two-gluon components we have 
calculated the flavor-singlet part of the 777 and 7/7 transi- 
tion form factors to NLO in a s and explicitly shown the 
cancellation of the collinear singularities present in the 
hard scattering amplitude with the UV one occurring in 
the unrenormalized distribution amplitudes. Assuming 
particle independence of the distribution amplitudes, we 
have employed the results for the transition form factors 
in an analysis of the available data 0, 0] and deter- 
mined the Gegenbauer coefficients to order n = 2 for 
the three remaining distribution amplitudes, the flavor 
octet, singlet and gluon one. The numerical results for 
the distribution amplitudes quoted for a = y/rif/CF are 
in agreement with the quark flavor mixing scheme pro- 
posed in [ill . 



The value for the lowest order gluonic Gegenbauer co- 
efficient is subject to a rather large error since the con- 
tributions from the two-gluon Fock components to the 
transition form factors are suppressed by a s as compared 
to the qq contributions. This suppression does not neces- 
sarily occur in other hard exclusive reactions; examples 
of such reactions, discussed by us briefly, are deeply vir- 
tual and wide-angle electroproduction of 7/ or 77' mesons 
as well as the (7*5*77(77') vertex. The latter two reactions, 
as it has turned out, are actually quite sensitive to the 
two-gluon components and future data for them should 
allow to pin down the gluon distribution amplitude more 
precisely than it is possible from the transition form fac- 
tor data. Other hard exclusive reactions which may be 
of relevance to our considerations are, for instance, the 
decays Xcj -> VV^'v' USUI or B —> rf'^K^ Last 
not least we would like to mention that the two-gluon 
components of other flavor-neutral mesons or even those 
of glueballs 01 can be studied in full analogy to the 77-77' 
case. 



VI. SUMMARY 



In this work we have investigated the two-gluon Fock 
components of the 77 and 77' mesons to leading-twist ac- 
curacy. Since the integral over the gluon distribution 
amplitude is zero, see l|2.4[) . there is no natural normal- 
ization of it in contrast to the case of the qq distribution 
amplitudes. Any choice of this normalization goes along 
with corresponding normalizations of the anomalous di- 
mensions and the projector of a two-gluon state onto a 
pseudoscalar meson. We have set up a consistent set of 
conventions for the three quantities which is imperative 
for leading-twist calculations of hard exclusive reactions 
involving 77 and/or 77' mesons. We have also compared 
this set with other conventions to be found in the litera- 
ture. 



Acknowledgments 



We wish to acknowledge discussions with M. Beneke, 
M. Dichl, T. Feldmann, D. Miiller and A. Parkhomenko. 
This work was supported by Deutsche Forschungs 
Gemeinschaft and partially supported by the Ministry of 
Science and Technology of the Republic of Croatia under 
Contract No. 0098002. 



15 



APPENDIX A: DEFINITIONS OF MESON 
STATES AND DISTRIBUTION AMPLITUDES 

The flavor content of the neutral pseudoscalar meson 
states we are interested in, is taken into account by 

7T° : -^=(uu-dd) -> C 3 = -^=A 3 , 

q% '■ -y=(uu + dd - 2ss) -> C 8 = -^=A 8 , (Al) 

qq 1 : —=(uu + dd + ss) — > Ci = If, 

V3 v^7 

where A^ are the usual 5f7 (3) Gell-Mann matrices and 1 
is the 3x3 unit matrix. For the flavor-singlet state, we 
use the general notation Q in which the flavor content 
is expressed in terms of n/ which denotes the number 
of flavors contained in qq 1 (nj = 3 in our case). This 



simplifies the comparison with the results for kernels to 
be found in the literature. 

As usual 0, 13 El we define the distribution ampli- 
tudes in a frame where the meson moves along the 3- 
direction. Neglecting the meson's mass its momentum 
reads 

p=\p+,0,0±], (A2) 

where we use light-cone coordinates v = [v + , v~ , vjj with 
v ± _ _|_ v 3y^/2 for any four-vector v 5 . We also 
introduce a light-like vector 

n = [Q, 1,0_l]. (A3) 

which defines the plus component of a vector, v + = n-v. 
The constituents of the meson, quarks or gluons, carry 
the fractions u and 1— it of the light-cone plus components 
of the meson's momentum. 



The distribution amplitudes are defined by Fourier transforms of hadronic matrix elements 



and 



fp 



2J2N, 



where z = [0, z , 0j_]. 



2ir 



*(-«) c t -^Ln^(z) 



'2N B 



P{p) 



fh 



:<t>Pg(u) 



- I" d>Z~~ j( u -(l-u))p-z n jL n v 



(n ■ p) J 2n 



G^{-z)SlG a »{z) P{p) 



(A4) 



(A5) 



Here, \I> denotes a quark field operator, G^ v the gluon 
field strength tensor, and its dual 



(A6) 



The quark and gluon operators in Eqs. (|A4|) . I|A5(1 are 

understood as color summed. The path-ordered factor 



exp < ig 



dsA(zs) ■ z 



(A7) 



where A is the gluon field, renders (j)pi and 4> p g g auge 
invariant. The distribution amplitudes in i|A4l IA5|) rep- 
resent either the unrenormalized ones {<Pp r ig {u)) if de- 
fined in terms of unrenormalized quark or gluonic com- 
posite operators or the renormalized one. In the lat- 
ter case the distribution amplitudes are scale depen- 
dent ((f>pi ig (u, /i 2 )). The distribution amplitudes defined 

Different conventions for the light-cone components are discussed 
in Ref. |47|. 



above satisfy the symmetry relations 

(j)p 1<s (u, \i 2 ) = 0pi,s(l - u, n 2 ) , 

(t) P g{u,H 2 ) = -<j> P g{\ ~ U,fl 2 ) . (A8) 

The definitions of the distribution amplitudes i|A4(l and 
(IA5(I can be inverted to 



in-pf P due^ 2u -^P- z Mu) , (A9) 
Jo 



and 



n^ny (0 



P 



G^ a (-z)nG a u (z) 
hn-pf V^/p I du e-^-Vf* <(>pg(u) 



(A10) 
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The projection of a collinear qq state onto a pseu- 
doscalar meson state is achieved by replacing the quark 
and antiquark spinors (normalized as vJ(p, X)u(p, A') = 
V2n-pS X y) by @ 



1 aj},rs,kl ~ u « 



>kl 



75 



(All) 



a/3 



where a (r, fc) and /3 (s, i) represent Dirac (flavor, 
color) labels of the quark and antiquark, respectively. 
When calculating amplitudes, the projector (|A11I) leads 
to traces. The projector holds for both incoming and 
outgoing states and corresponds to the definition of the 
the quark distribution amplitudes HA4|> . It is to be used 
in calculations of hard-scattering amplitudes which are 
to be convoluted with fp/{2y/2N c )4>pi subsequently. 

The form of the projection of a gg state on a pseu- 
doscalar state with momentum p can be deduced by not- 
ing that the helicity zero combination of transversal gluon 
polarization vectors can be written as 0] 

e*(up, A)e"((l - u)p, -A) - e^(up, -A)e"((l - u)p, A) 
= i sign(A) , 

(A12) 

where e 12 = —e 21 = 1 while all other components of 
the transverse polarization tensor are zero. It can be 
expressed by 



n ■ p 



(A13) 



Instead of n any other four vector can be used in (|A13|) 
that has a non-zero minus and a vanishing transverse 
component. The projector of an state of two incoming 
collinear gluons of color a and b and Lorentz indices \i 
and v, associated with the momentum fractions u and 
(1 — u), respectively, onto a pseudoscalar meson state 
reads 



V 9 h = 



2 y/N* - 1 u(l - u) 



(A14) 



The complex conjugated expression is to be taken for an 
outcoming gg state. The projector is to be used along 
with the distribution amplitude fp/(2y / 2N c )(j>p g . The 
additional factor [u(l — m)] _1 appearing as part of the 
projector, is a consequence of the fact that in pertur- 
bative calculations of reactions involving two-gluon Fock 
components, the potential A of the gluon field occurs, 
while the gluon distribution amplitude is defined in terms 
of the gluon field strength operator, see l|A5|l . The con- 
version from a matrix element of field strength tensors 
(|AT0|) to one of potentials is given by |32j, l49j 



(0\A a {-z)A ,3 (z)\P) 




du e" 



i(2u-i)p-z ^Pgiu) 
u(l — u) 



(A15) 



The gluonic projector i|A14(l is obtained (up to the fac- 
tor [u(l — it)]^ 1 explained above) by the coupling of two 



collinear gluons into a colorless pseudoscalar state. In the 
context of mixing under evolution another normalization 
of it appears to be more appropriate, see l|2.41|l . This 
normalization is accompanied by corresponding changes 
in the gluon distribution amplitude <pp g and the anoma- 
lous dimensions, as is discussed in detail in Sec. II. 
For Levi-Civita tensor we use the convention 

e 0123 = -1 , (A16) 

which leads to 

Tr [757 / V7 t V] = Aie^ 13 (A17) 
(with 75 = «7 7 1 7 2 7 3 ). 



APPENDIX B: THE P-y TRANSITION FORM 
FACTOR - 
DETAILS OF THE CALCULATION 

In this appendix, we provide some details of the cal- 
culation of the evolution kernels and the hard scattering 
amplitude for the flavor-singlet contribution to the P7 
transition form factor. These quantities can, in princi- 
ple, be taken from the literature (see, e.g. 0, Q and 
|50| 6 ) but the conventions and notations differ. How- 
ever, since it is imperative to use a consistent set of con- 
ventions for the hard scattering amplitude and the dis- 
tribution amplitudes, we recalculate them. In doing so 
we follow closely Ref. 01 ■ Dimensional regularization 
in D = 4 — 2e dimensions is used to regularize UV and 
collinear singularities which appear when calculating the 
one-loop diagrams. According to 0, the 75 problem, 
i.e, the ambiguity which enters the calculation due to the 
presence of one 75 matrix and the use of dimensional reg- 
ularization method, is resolved by matching the results 
for the hard-scattering part with the results for the per- 
turbatively calculable part of the distribution amplitude, 
since the physical form factor is free of ambiguity. We 
employ the MS coupling constant renormalization along 
the same lines as in [lj]. We note in passing, that as 
long as the singularities are not fully removed from the 
amplitudes, the following relations are to be used for the 
change of the scale of the coupling constant 



a s (/u 2 



M 2 



,(H 2 R )[l + 0(a s )] (Bl) 



and for the (3 function 
f3(a s (p 2 ), e) = fi 2 -^a s (fi 2 ) 



-ea s (/i 2 ) 



« 2 (M 2 ) 

47T 



A)- 

(B2) 



In Ref. |5Cl the NLO corrections to the deeply virtual Compton 
amplitude ■y'p —* -y*p have been calculated. In the limiting case 
of zero skewness the Compton amplitude is related to our process 
by crossing. 
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FIG. 12: Sample NLO Feynman diagrams contributing to the 7*7 — > qq amplitude. 




Gl 



G2 



fCRfiJWiRP 



G3 



FIG. 13: Distinct one-loop Feynman diagrams contributing to the 7*7 — + gg amplitude. Other contributing diagrams are 
obtained from these by reversing the direction of the fermion flow in the loops. 



The usual renormalization group coefficient is given by 



0o 



11 



-, n f 



(B3) 



1. Amplitudes 



The amplitude 77 — * qq denoted by T q g (examples of 
contributing Feynman diagrams are depicted in Fig. I12|) 
has the structure already quoted in (|2.15|l where 



1 



1 

u 



T$(u) = —Ail\Ju)+A^(u). (B4) by 



-1 



,(i) 



col,qq 



(D, 



The functions A read 
1 



1 



3 + 21n(l -u) + 1 -u), 



1 



1 - 7 
+0 



1 -it 



ln(l - u) + ln 2 (l - u) 



1 - u) . 



(B5) 



In obtaining the above results the projector (|A11() is em- 
ployed. The results for the flavor-octet and singlet cases 
differ only in the flavor factors (see (|2.17|l and 1|2.44J0 . 

Next, we calculate the amplitude T gg for the subpro- 
cess 7*7 — > gg. The appropriate gluonic projector is the 
complex conjugate of l|A14() . For the case of the transi- 
tion form factor we can work in a Breit frame where the 
momentum of the real photon, q2, is proportional to the 
vector n from Eq. 1|A3() . and can therefore be employed 



in i|A13l) . There are 6 one- loop diagrams that contribute 
to this subprocess amplitude. Three representative di- 
agrams (Gl, G2, G3) are shown in Fig. ^| The other 
three reduce to the first three ones by reversing the di- 
rection of the fermion flow in the loop. Moreover, it is 
easy to see that 



T G2 = -T G1 (u -> 1 - u) . 



(B6) 



Thus, one has only to calculate the contributions from 
the diagrams Gl and G3. 

The complete unrenormalized NLO contribution is the 
sum of individual contributions in which, expectedly, the 
UV singularities cancel. The hard-scattering amplitude 
Tgg has the structure quoted in l|2.15|l where T, 



(i) 

g-g' is given 



gg 



-A 



(i) 

coLgg 



(B7) 



The functions A read 7 



A 



(i) 

col,gg 



(«) 



1 



(«) = 



1t(l — It) 

1-71 



ln(l — u) — (u — ► 1 — u) 



3 - - ln(l - u) 



2 tt 



ln 2 (l -u)-(u->l-u) 



(B8) 



7 Making use of the crossing relations, it can be shown that the 
functions IB5I and IB8I are in agreement wit h th e coefficient 
functions for the Compton amplitude quoted in 15(1 . 
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Dl 



D2 



FIG. 14: LO Feynman diagrams that contribute 
to <j> qg . The crossed circle denotes the vertex of 
{0\^(-z)C 1 fr f5 /V2K' i f(z). 



2. Kernels 

For the calculation of the renormalization matrix Z, 
respective we utilize the method proposed 

in [12l EE) of saturating the mcsonic state by its valence 
Fock components (I2.1|l which leads to 



$p r (u) = -i<t>(u,v) <g) 




(B9) 



The elements of the matrix 4> are defined as in l|A4|l and 
(|A5|I with the replacement of |P) by \qq±) and Igg). They 
are thus perturbatively calculable and determine the ma- 
trix Z . 

The calculation of the matrix element Z qq proceeds 
along the same lines as indicated for the flavor-octet case 
in Ref. and the contributing diagrams are displayed 
there. The respective kernel V qq reads 



V qq (u,v) = 2C f 



1 



1 



@(v-u) 




(BIO) 



where the usual plus distribution is defined as 



[f(u,v)} =F(u,v)-S(u-v) j dzF{z,v). (Bll) 

This result also holds for the flavor-octet case. 

We proceed to the evaluation of Z qg , or rather V qg . 
According to the definition of the qq 1 distribution ampli- 
tude, the matrix element that is of interest here, is given 
by (z=[0,z-,O±]) 



<t>qg{u) = 



dz 
~2tt 

( 



i(2u-l)p-; 



fill 



(B12) 



The relevant Feynman diagrams for the calculation of 4> qg 
are depicted in Fig. The qq vertex, (g), is of the form 



Ci- 



S(u n ■ p — n ■ k) 



(B13) 



where k represents the momentum of the quark entering 
the circle. The vertex (|B13() occurs also in the calcula- 
tion of the 4> qq where the LO contribution is obtained by 
contracting the vertex just with the qq projector (|A11() 
and, hence, one obtains (f> qq (u, v) — 5(u — v) as it should 
be (see l|27T3)ll. 

Due to the presence of only one 75 matrix, we are 
confronted with the 75 problem, as in the calculation 
of T qq . When using the naive 75 scheme, in which the 
75 matrix retains its anticommuting properties in D 
dimensions, we obtain three different results depending 
on the position of 75 inside the trace: 



(j>qg,Dl(u,v) 



(4tt) 



2c 



d D l 



4tt 1 1 



(l-vf 



2c 



11 



1 - e/2" \v 2 (l - v) 
(1-u) 



&{v - u) 



v(l — v) : 



T Q(u-v) 



where 



*e{-(2«-i),-i,i}. 



(B14) 



(B15) 



The loop integral can be worked out analytically 8 and we 
refer to for the result. 
One can easily see that 

4>qg,D2 ( U i V ) = -4>qg,Dl (u, 1 ~ v) , (B16) 

and finally 

<t>qg(u,v) = 4>qg,Dl(u,v) - <j) qg .Dl(u, 1 - v) . (B17) 

The kernel V qg is a residue of the UV singularity embod- 
ied in the loop integral appearing in (|B14J| and, hence, 
is related to the term multiplying the integral in (|B 14|) . 
Since the term proportional to S is finite (~ e(l/e)), it 
does not contribute to V qg . Moreover, since (j) qg being 
antisymmetric under the replacement of v by 1 — v, is 
to be convoluted with the matrix element (gg\P) (see 
(|B9|n . which has the same symmetry properties as the 



Th e tr eatment of the integral in Eq. IB14I was explained in detail 
in |l2j. The crucial point is to retain a distinction between UV 
and collinear singularities. 
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full gluon distribution amplitude (see one can re- 

place <j> qg by (j/ (u, v) = 2(f>qg t Di(u, v) in order to obtain 
a more compact representation of the kernel 



[u(l — u)] 1 is included in <f>p g instead in the gg projec- 
tor (JSEJl. The results for T gg l|H7|. (Eg) will, hence, be 
multiplied by u(l — u), while the kernels take the form 



V qg (u, v) = -2^Jn f C F < -^@( v -u)- 



u 
v 



(B18) 



The set of LO evolution kernels is completed by 



V gq (u,v) = 2y/n f C F i —Q(v-u)- 



Vgg(u,v) = 2N C 



1-u 
1-u 



(B19) 



f« 


[( 


Q(v 


I" 




v - 


u — > 


1 


— u 


u — > 


1 


— V 



2u- 1 



Q(v-u) 



(3 5(u - v) 



(B20) 



Since, except of the normalization, there is general agree- 
ment in the literature on these kernels, see e.g. 0,U, we 
quote them without giving any detail of their calculation. 
Finally, we comment on an alternative definition of the 
gluon distribution amplitude which one occasionally en- 
counters in the literature. In that definition the factor 



V. 



qg V qg v(l - v) , 



v„, 



Va, 



V a , 



„(!-„) U(1 ~ U) W 
u(l — u) 



The result for the transition form factor, as for any other 
physical quantity, is, obviously, invariant under the re- 
definition of the gluon distribution amplitude. 



APPENDIX C: SOME PROPERTIES OF THE 
EVOLUTION KERNEL 



It is easy to verify that the evolution kernels (|B10|) and 
(|B18|) - (IB20() satisfy the symmetry relations 



V (I-V) V qq {u,v) 
V 2 {l- V f V gg (u,v) 



U (1-u) Vqq{v,U) , 

U 2 (l - Uf Vgg(v,u) , 



\l-v) 2 V qg (u,v) = u (1-u) V gq (v,u). (CI) 



J 



The kernels Vij, convoluted with the weighted Gegenbauer polynomials C™ of order m = 3/2, 5/2, result in 



Vqq(u 


v) C 


5 v (1 - 


v) C*J 2 {2v- 


1) 


— ^11 

In 


U (1- 


u) C^ 2 {2v 


-1), 


V qg (u, 


v) (5 


5 v 2 {l- 


v) 2 C 5 S\(2v- 
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— ^19 
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U (1- 


«) cT(2v 


-1) 


V gq (u 


v) i 


5 v (1 - 


v) C?/ 2 (2v- 
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— n aq 

In 


u 2 (l - 


u?C h J\{2v 


-1), 


Vgg{ U , 


v ) <i 


5 v 2 {\- 


v) 2 C 5 S\(2v - 


1) 


— ^,99 
In 


u 2 (l- 


u) 2 Ct\(2v 


-1) 



(C2) 



The factors on the right hand side of (|C2|) multiplying the 
Gegenbauer polynomials are the anomalous dimensions. 
The results quoted for them in (|2.28|l can be read off 
from (|C2J) (for a detailed discussion see Q). Finally, we 
mention that the off-diagonal anomalous dimensions in 
(|2.28() satisfy the relation 



In 
In 



N 



5/2 



3/2 



where 



3/2 (n + l)(n + 2) 
" 4(2n + 3) ' 

M 5/2 »(n + 3) 3/a 
36 " ' 



(C3) 

(C4) 
(C5) 



represent the normalization constants of the correspond- 
ing Gegenbauer polynomials 



duu{l-u)C^ 2 {2u-l)Cl' 2 {2u-l) = N^ 2 5 nm , 

duu 2 (l~u) 2 C* /2 {2u-l)CU 2 {2u~l) = N^ 2 S nm . 

(C6) 



Throughout the paper we investigate only the LO behav- 
ior of the evolution kernels and corresponding anomalous 
dimensions. Beyond leading order, the relations corre- 
sponding to 1|C2|1 and (|C3I) get modified due to mixing of 
conformal operators starting at NLO (see, for example, 
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